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Abstract. We study dispersive properties for the wave equation in the Kerr space- 
time with small angular momentum. The main result of this paper is to establish 
uniform energy bounds and local energy decay for such backgrounds. This follows a 
similar result for the Schwarzschild space-time obtained in earlier work |21| by the 
authors and collaborators. 



1. Introduction 



The aim of this article is to study the decay properties of solutions to the wave equation 
in the Kerr space-time, which describes a rotating black hole. Until recently even the 
problem of obtaining uniform bounds for such solutions was completely open, and only 
some partial results were obtained in [TJ . We also refer the reader to related independent 
work in [13]. However, the techniques used in these papers are of a different flavor, as 
they do not carry out such a precise analysis of the dynamics near the trapped set. 

Our aim here is to establish global in time energy bounds for the wave equation in the 
Kerr space-time, as well as a local energy decay estimate. These bounds apply in the full 
region outside the event horizon, as well as in a small neighborhood on the inside of the 
event horizon. 

The starting point in our analysis is the earlier work [21] of the authors and collabora- 
tors, which establishes similar bounds for the wave equation in the Schwarzschild space- 
time. The idea is to treat the Kerr geometry as a small perturbation of the Schwarzschild 
geometry, and then adapt the methods in [21]. Consequently in this article we are only 
considering Kerr black hole backgrounds with small angular momentum, which are close 
to the Schwarzschild space-time. Nevertheless, we are confident that our methods will 
carry over also to the case of large angular momentum. 

Another goal of the earlier article [21] was to establish Strichartz estimates in the 
Schwarzschild space-time. We also aim to consider the similar problem for the Kerr 
space-time. However, this requires very different technical tools, and will be considered 
in a subsequent paper. 

The local energy estimate in [21] is proved using the multiplier method; the delicate 
issue there is to show that a suitable multiplier can be found. This method is quite robust 
under small perturbations of the metric, and for the most part it easily carries over to the 
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Kerr backgrounds with small angular momentum. There is however one region where this 
does not apply, precisely near the photon sphere r = 3M (which contains alQ the trapped 
periodic geodesies in the Schwarzschild space-time) . Hence most of the new analysis here 
is devoted to understanding what happens there. 

The paper is organized as follows. In the next section we discuss the classical lo- 
cal energy decay decay estimates in the Minkowski space-time and small perturbation 
thereof. Then we provide a brief overview of the local energy estimates proved in [21] 
for the Schwarzschild space-time, along with a discussion of the relevant geometrical is- 
sues. Finally, the last section contains the description of the Kerr space-time and all 
the new results. Our main local energy estimate is contained in Theorem 14.11 This is 
complemented by higher order bounds in Theorem 14.41 



2. Local energy decay in the Minkowski space-time 

In the Minkowski space-time R^^^, consider the wave equation with constant coeffi- 
cients 

(2.1) Du = /, m(0) = Mo, dfu{0) = ui 
Here □ = — A. More generally, let 

be the usual d'Alembertian associated to a Lorentzian metric g. 

The seminal estimate of Morawetz asserts that for solutions to the homogeneous 
equation Du = we have the estimate 

(2.2) / / -^\fuWt,x)dxdt+[\u{t,0)\''dt<\\Vuo\\l2 + \\ui\\l2 

JR JR'^ Fl JR 

where y/ denotes the angular derivative. This is obtained combining energy estimates 
with the multiplier method. The radial multipHer Qu = {dr + is used, where r 
denotes the radial variable. 

Within dyadic spatial regions one can control the full gradient Vu, but the square 
summability wth respect to dyadic scales is lost. Precisely, we define the local energy 
norm 

(2.3) \\u\\lEm = sup2"^||u||i2(Rx{|^|e[2J-i,2J]}) 
and its counterpart 

(2.4) MleI = \\^x,tu\\LE„ + \\\x\'^u\\lem 
For the inhomogeneous term we use the dual norm 

(2.5) ll/IU-EJf = 22 ||/||^2(Rx||:t|g[2'°-l,2''])) 

fcez 



except of course for the rays along the event horizon, which are not relevant to this discussion 
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Then we have the following scale invariant local energy estimate for solutions u to the 
inhomogeneous equation l|2.ip : 

(2.6) llV^ilUrL^, + \MleI, < llVliolU^ + \\ui\\l2 + \\.f\\LEl,+LlLl 

This is proved using a small variation of Morawetz's method, with multipliers of the form 
a{r)dr + b{r) where a is positive, bounded and increasing. 

There are many similar results obtained in the case of perturbations of the Minkowski 
space-time; see, for example, [19], [H], [16], [23], [21], [25], [2], [2Q]. Relevant to us is the 
case of small long range perturbations of the Minkowski space-time, considered in |22j . 
The metrics g in considered there satisfy 

(2.7) V sup \g{t,x)^gM\ + \x\\VxMt,x)\ + \x\^\Vltg{t,x)\<e 

where gM stands for the Minkowski metric. Then as a special case of the results in [22] 
we have 

Theorem 2.1. [22j Let g be a Lorenzian metric in M'^+^ which satisfies (|2.7p with e small 
enough. Then the solution u to the inhomogeneous problem 

(2.8) = u(0) = uo, dtu{Q) = ui 
satisfies the estimate (|2.7p . 

No general such results are known for large perturbations, where on one hand trap- 
ping for large frequencies and on the other hand eigenvalues and resonances for small 
frequencies create major difficulties. The Schwarzschild and Kerr metrics are such large 
perturbations where trapping plays a major role. 

3. Local energy decay in the Schwarzschild space-time 



In the original coordinates the Schwarzschild space-time is given as a metric whose 
line element is (for / = M x (2M, oo) x S^) 

(3.1) ds^ ^-{^- V")"^^' + " V") + 

where dw^ is the measure on the sphere and t, r are the time, respectively the radius 
of the §^ spheres. This metric is well defined in two regions, 

/ = R X (2M, oo) X §^ // = M X (0, 2M) x 

Let Us denote the associated d'Alembertian. 

The singularity at r = is a true metric singularity. However, the singularity at the 
event horizon r — 2M is an apparent singularity that can be removed by a different choice 
of coordinates. Following [15], let 

r* = r + 2M log(r - 2M) - 3M - 2M log M 

and lei V = t-\-r* . In the new coordinates (r, v, u) the metric becomes 

ds^ = - [l- — ^du^ + 2dvdr + r'^dto'^ 
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and can be extended to a larger manifold / U //. Moreover, if w = t — r* , one can 
introduce global nonsingular coordinates (all the way to r = 0) by rewriting the metric 
in the Kruskal-Szekeres coordinate system, 

v' = e™", w' = —e^™. 

There are two places where trapping occurs on the Schwarzschild manifold. The first 
is at the event horizon r = 2M, where the trapped geodesies are the vertical ones in the 
(r, v.ijj) coordinates. However, this family of trapped rays turns out to cause no difficulty 
in the decay estimates since the energy decays exponentially along it as w ^ oo. The 
second family of trapped rays occurs on the surface r = 3M which is called the photon 
sphere. Null geodesies which are initially tangent to the photon sphere will remain on the 
surface for all times. Unlike the previous case, the energy is conserved for waves locaHzed 
along such rays. However, what makes local energy decay estimates at all possible is the 
fact that the trapped rays on the photon sphere are hyperboHc. 

The {r,v,ijj) coordinates are nonsingular on the event horizon, but have the disad- 
vantage that the level sets of v are null surfaces. This is why it is more convenient to 
introduce 

V = V — 

where /i is a smooth function of r. In the (-u, r, w) coordinates the metric has the form 

ds^ = - (l - + 2 (^1 - (l - dvdr 

+ (2/i'(r) - (l - '^Y^,'{r)fyr'+r^duj\ 

On the function /i we impose the following two conditions: 

(i) ^(r) > r* for r > 2M , with equality for r > 5M /2. 

(ii) The surfaces v = const are space-Hke, i.e. 

M'(r)>0, 2- (l-^)Ai'(r) >0. 

The first condition (i) insures that the {r,v,uj) coordinates coincide with the {r,t,uj) 
coordinates in r > 5M/2. This is convenient but not required for any of our results. 
What is important is that in these coordinates the metric is asymptotically fiat as r ^ oo 
according to l|2.7p . 

Given < re < 2M we consider the wave equation 

(3.2) Osu = / 
in the cylindrical region 

(3.3) MR^{i>0, r> rj 
with initial data on the space-like surface 

(3.4) = Mr n {i = 0} 
The lateral boundary of Mr, 

(3.5) J:+^MRn{r = r,} 



LOCAL ENERGY ESTIMATE ON KERR BLACK HOLE BACKGROUNDS 



5 



is also space-like, and can be thought of as the exit surface for all waves which cross the 
event horizon. 

We define the initial (incoming) energy on T,~^ as 



(3.6) 



{\dru\^ + \di,u\^ + \fu\^) r^drduj 



the outgoing energy on as 



(3.7) 



and the energy on an arbitrary v slice as 



(3., 



E[u]{io) = 



{\dr 



\dr,uY + \yiuY)ridMuj 



The choice of the local energy norm LEs is inspired from (|2.3p . However, there is 
a loss along the trapped geodesies on the photon sphere. Consequently, we introduce a 
m odifiecil 

local energy space 



(3.9) 



M\les 



LE^ 



and local energy space 

(3.10) IIwIIlbi = \\drU\\LEM + \\dvU\\LEs + \\fu\\LEs 

For the inhomogeneous term we use the norm 



\r U\\LE^ 



(3.11) 



\LE% 



1 - 



3M 



In the three formulas above we implicitely assume that all norms are restricted to the set 
Mr where we study the wave equation l|3.2p . Then we have the following result: 



l/ll 



Theorem 3.1. [21] Let u he so that Dsu = /. Then we have 



(3.12) 



£;[u](S]+) + sup£:M(s) 



\LE 



l.<E[u]{-^-j 



Note that, compared to the norms LEm-, LEIj, the weights have an additional poly- 
nomial singularity at r = 3M, but there are no additional losses at the event horizon or 
near oo. Furthermore, by more refined results in [21], this polynomial loss can be relaxed 
to a logarithmic loss, i.e. the factor 1 — can be improved to | ln(r — 3Af)|~^ near 
r = 3M. This is related to the fact that the (periodic) trapped rays on the photon sphere 
are hyperbolic. 

We also remark that in the expression of LEg it was sufficient to measure drU. This 
is due to the implicit cancelation caused by the fact that the symbol of the operator dr 
vanishes on the trapped set. 



^notations are slightly changed compared to |21| in order to insure some uniformity across the three 
models described in the present paper 
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The choice of £ (0, 2M) is unimportant since the r-shces r = const e (0, 2M) 
are spaceHke. Hence moving from one such r-shce to another is equivalent to solving a 
local hyperbolic problem, and involve no global considerations. Thus in the proof of the 
theorem one can assume without any restriction in generality that re is close to 2M . 

Local energy estimates were first proved in [18] for radially symmetric Schrodinger 
equations on Schwarzschild backgrounds. In [Hill [5], those estimates are extended to 
allow for general data for the wave equation. The same authors, in [61I2, have provided 
studies that give certain improved estimates near the photon sphere r = 3M. Moreover, 
we note that variants of these bounds have played an important role in the works [8] and 
[TT] . [12] which prove analogues of the Morawetz conformal estimates on Schwarzschild 
backgrounds. 



4. Local energy decay in the Kerr space-time 



The Kerr geometry in Boyer-Lindquist coordinates is given by 

2 , „ ,.11 , „ J„2 , „ . ,2 , „_jq2 



ds = gttdt + gt^dtdcj) + grrdr + 

where t G R, r > 0, (0, 9) are the spherical coordinates on §^ and 

A - sin^ 9 2Mr sin^ 9 
9tt = 5 , gt<p = -2a , grr = -r- 

" - sm 9, ggg = p 



p2 



with 



A = - 2A/r + a^, p^ = + cos^ 6*. 

A straightforward computation gives us the inverse of the metric: 

_ (r2 + a^f ~ a^Asin^ 9 _ 2Mr „. _ A 

" ^ ^"a ' ^ ~ ~" ^2 A ' ^ ^72- 



p2A 



A - sin^ 9 _ 1 

p2Asin2 6' ' ^ ^ p2 



The case a = corresponds to the Schwarzschild space-time. We shall subsequently 
assume that a is small a <C Af, so that the Kerr metric is a small perturbation of the 
Schwarzschild metric. We let denote the d'Alembertian associated to the Kerr metric. 

In the above coordinates the Kerr metric has singularities at r = on the equator 
9 = it/2 and at the roots of A, namely r± = M ± V — a^. As in the case of the 
Schwarzschild space, the singularity at r r_|_ is just a coordinate singularity, and corre- 
sponds to the event horizon. The singularity at r = r_ is also a coordinate singularity; 
for a further discussion of its nature, which is not relevant for our results, we refer the 
reader to [lOl [15] . To remove the singularities at r = r± we introduce functions r* , v+ 
and so that (see [15]) 

dr* ^ {r^ + a^)/\~'^dr, dv+ = dt + dr*, d(t>+ = dcf) + aA~'^dr. 
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The metric then becomes 



2 Air 

ds^ = - (1 ^)dv'i + 2drdv+ - ^^ap-'^Mrsir? edv+d(j)+ - 2asin2 9drd(j)+ + p^dO'^ 



+ + a^f - Aa^ sin^ 0] sm^ ^-^-^^ 



which is smooth and nondegenerate across the event horizon up to but not including 
r = 0. Just like in [21], we introduce the function 

V ^ v+ — p{r) 

where is a smooth function of r. In the (w, r, (j)+,9) coordinates the metric has the form 

ds' - (1 - ^)dv' + 2 (^1 - (1 - ^)fi'{r)] didr 

2/i'(r)-(l —){^'{r)fyr^ 



- 2a0{l + 2p-'^Mrn'{r)) sin^ drd(t)+ + p^dO'^ 
+ p"^[(r2 + - Aa^ sin^ 9] sin^ 6ld0^. 

On the function /i we impose the following two conditions: 

(i) p{r) > r* for r > 2M, with equality for r > 5A//2. 

(ii) The surfaces v = const are space-like, i.e. 

2Mr 

p'ir)>0, 2-{l-^)p'{r)>0. 

P 

As long as a is small, we can work with the same function p as in the case of the 
Schwarzschild space-time. 

For convenience we also introduce 

CW<^+ + (1 - CW)'/' 

where C is a cutoff function supported near the event horizon and work in the (u, r, </>, 9) 
coordinates which are identical to {t, r, <f>, 9) outside of a small neighborhood of the event 
horizon. 

Carter [9] showed that the Hamiltonian flow is completely integrable by flnding a 
fourth constant of motion K that is preserved along geodesies. If E and L are the two 
constants of motion associated with the Killing vector flelds dt and 9^, the equations for 
the null geodesies can be reduced to the following (see, for example, [lO] or [26] ) 

p t ~ a[L — £/asm 9) H ^ 

L- Easin^e {r"^ + a'^)aE - L 



4^2 {L~Eas\T?ef 



(4.1) sin^e* A 



p-'9' ^K- 2 

sin 



= - /^A + ((r^ + a^)E - ahf 

where the overdot denotes differentiation with respect to an afHne parameter s. This 
parametrization of the null geodesies is nondegenerate away from the surfaces r = r±. 
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Next we discuss the geometry of the trapped null geodesies. The level sets r = tq of 
r are time- like for tq > null for r = r+ and space-like for r_ < tq < r+. The latter 
implies that there are no trapped null geodesies inside the region {r_ < r < r+}. On the 
null surfaces r = r±, through each point there is a unique null vector which is tangent 
and which generates a trapped null geodesies. 

To find the trapped null geodesies in the region r > it suffices to consider the 
behavior of the fourth degree polynomial 

P(r) = -A'A + ((r^ + a^)E - alf 

in the last equation in l|4.ip . At least one of the parameters E, K and L should be 
nonzero, and the third equation shows that K > Q and that we cannot simultaneously 
have E = K = 0. Thus P is always nondegenerate. The key observation is that that the 
simple zeroes of P correspond to turning points in the last equation, and only the double 
zeroes are steady states. There are several cases to consider. 

a) If = then K > 0. Thus P has at most one positive root, where it changes sign 
from + to — . This root is a right turning point for the ode, and there are no trapped 
null geodesies. 

h) E ^ 0. Then P has degree 4 and P > in [r_,r+]. If P has any zero in [r_,r+] 
then the square expression must vanish, and this zero must be a double zero. We claim 
that in (r+,oo) P has either no root or two roots (counted with multiplicity); this is 
easily seen, as P must have either (at least) two complex conjugate roots or a negative 
root (the sum of the roots equals 0) and (at least) another one smaller than r_ (since 
P(r_) > 0). There are three subcases: 

bl) P has no roots larger than r+. Then r is monotone along null geodesies, and there 
are no trapped null geodesies. 

b2) P has two distinct positive roots r+ < ri < r2 . There it must change sign from + 
to — , respectively from — to +. Hence ri is a right turning point and 7'2 is a left turning 
point for the ode. Thus no trapped null geodesies exist. 

b3) P has a double positive real root tq. Then this root is a steady state, and all other 
solutions converge to the steady state at one end, and escape to or infinity at the other 
end. 

This analysis shows that the only trapped null geodesies are those along which r is 
constant. The polynomial P has a double root if the following two relations hold, 

((r^ + a^)E - aLf = KA, 2rE{{r^ + a^)E - aL) = K{r - M) 

which we rewrite in the form 

r^E^A ^ / n o 2rA \ 

The right hand side in the 9 equation must be nonnegative. Substituting in the above two 
relations we obtain a necessary condition for the existence of trapped geodesies, namely 
the inequality 



(4.2) 



(2rA - (r - M)p2)2 < 4aV2 A sin^ 9 
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One can show that this condition is also sufficient. The expression on the left has the 
form 

2rA - (r - M)p^ =r^{r- 3M) + 2ra^ - {r - M)a^ cos^ 6 

If a = then it has a single positive nondegenerate zero at r = 3M, which is the photon 
sphere in the Schwarzschild metric. Hence if < a ^ Af it will still have a single zero 
which is close to 3M. A rough computation leads to a bound of the form 

(4.3) |r-3A/|<2a, a < 2Af 

Thus all trapped null geodesies lie within 0{a) of the r = 3M sphere. 

We would Hke a characterization of the aforementioned trapped geodesies in the phase 
space. Let t, ^, $ and 6 be the Fourier variables corresponding to t, r, and 6, and 

p{r, 0, r, e, $, e) = .g"r2 + 2g*^T^ + g**^^ + g'^e + 9^'^^ 

be the principal symbol of Dk- On any null geodesic one has 

(4.4) p(i,r,(/),0,r,^,$,e) = O. 



Moreover, the Hamilton flow equations give us 

dp _ 2A _ 



(4-5) r^-^^-^i 



(4.6) i = ^^= g^y + 2g'fr<^ + g''f<^^ + g^^f + g'^Q' 
We rewrite the latter in the form 

(4.7) pH - ^ „2i?,(r, r, $)A-2 + p'd,{p-^)p + 2(r - M)f 
where 

i?a(r,T,$) = {r^ + a^){r^ -~ Mlr^ + a^r + M)t^ - 2aM{r^ - a^)r$ - a^(r - A/)$^ 



For geodesies with constant r, one needs to impose the additional condition r = 0. 
Hence from l|4.5p either r = r±, which correspond to the geodesies at r = 2A'/ in the 
Schwarzschild case, or ^ = 0. In the latter case from (|4.7p we obtain a polynomial 
equation for r, namely 

(4.8) Ra{r,T,^) = 

Furthermore, due to (|4.4p we must also have the inequality 

-((r2 + a2)2 _ a^Asin^ 6^ - 2aA/rT$ + ^ ~ < 

sin 9 

If a is small and r is as in (|4.3p this allows us to bound $ in terms of r, 

(4.9) |$| < 4Af|r| 

For <f> in this range and small a the polynomial T^^Ra{r, r, $) can be viewed as a small 
perturbation of 

r"2i?o(r,r,$) = r^(r - 3Af) 
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which has a simple root at r ~ 3M. Hence for small a the polynomial Ra has a simple 
root close to 3M, which we denote by ra{T,4>). By homogeneity considerations and the 
implicit function theorem we can further express ra in the form 



r,(T,$) = 3A/f , re C°°([-e,6] x [-4,4]) 

Since ro(r, $) = 3M, it follows that we can write ra{T, $) in the form 

ra{r,<^) = 3M + aF(^-^,^y F G C-([-e, e] x [-4, 4]) 

The above analysis shows that the trapped null geodesies corresponding to frequencies 
(t, $) are located at radial frequency ^ = and position r = rQ(r, ^). One would be 
naively led to define the local smoothing spaces associated to the Kerr space-time by 
replacing the factor r — 3M in l|3.9p and l|3.1ip with the modified factor r — ra(T, $). 
Unfortunately, this is no longer a scalar function, but a symbol of a pseudodifferential 
operator. In addition, this operator depends on the time Fourier variable r, which is 
inconvenient for energy estimates on time (v) slabs. 

Consequently, we replace the r — ra {T, $) weight with a polynomial in r which has the 
same symbol on the characteristic set p = 0. More precisely, for r close to 3M we factor 

p(r, 0, T, e, $, e) = .g"(T - Ti(r, 0, C, $, e))(T - T2(r, 0, ^ $, 9)) 

where ri, T2 are real distinct smooth 1-homogeneous symbols. On the cone t — Ti the 
symbol r ~ ra (r, (/)) equals 

c,(r, 0, e, $, e) = r - ra(r„ $) = r - MI - aF (^^, , i = 1, 2 

If r is close to 3M and |a| ^ Af then on the characteristic set of p we have \4>\ < 4A'/|r|, 
therefore the symbols Ci are well defined, smooth and homogeneous. They are also 
nonzero outside an 0{a) neighborhood of 3M. 

We use the symbols q to define associated microlocally weighted function spaces L^. in 
a neighborhood / x §^ of 3A/ x which does not depend on a for small a. For functions 
u supported in / x §^ we set 

2 



There is an ambiguity in this notation as we have not specified the coordinate frame 
in which we view Ci as a pseudodifferential operator. However, it is easy to see that 
different frames lead to equivalent norms. We note that the low frequencies in a are also 
irrelevant, and can be removed with a suitable cutoff. After removing the low frequencies, 
the quantization that we use for Ci becomes unimportant as well. We also define a dual 
norm aL^ for functions g supported in / x namely 

MIl^^ , inf (Il.9i|li2 + 11.92111,0 

Ci{x,D)gi+g2=g 



Since the symbols q are nonzero outside an 0{a) neighborhood of 3Af , it follows that 
both norms L^. and are equivalent to outside a similar neighborhood. 
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Now we can define local smoothing norms associated to the Kerr space-time. Let x(r) 
be a smooth cutoff function which is supported in the above neighborhood / of 3M and 
which equals 1 near 3M. Then we set 

(4.10) 

Mle], =\\xiDt - T2{D,x))xu\\l2^ + \\x{Dt ~ n{D,x))xu\\Li^ + ~ X^)dtu\\LEM 

+ 11(1 - X^)fu\\LEM + \\dM\LE^, + \\r^^u\\LEM 

We remark that this norm is equivalent to the Minkowski norm LE\j outside an 0{a) 
neighborhood of 3M, but it is degenerate on the trapped set. 

For the nonhomogeneous term in the equation we define a dual structure, 

\\f\\LE*, 11(1 - X)/I|lSJ, + \\xf\\c^L^+C2L^ 

To state the main result of this paper we use the notations in (|3.3p - p.8p . with the 
parameter chosen so that r_ < r^, < 

Theorem 4.1. Let u solve = / in Mb- Then 

(4.11) ll^^lli^, + supi?M(«) + E[u]{Y.+) < i?N(S;,) + ll/llis.,. 

V 

in the sense that the left hand side is finite and the inequality holds whenever the right 
hand side is finite. 

The proof of the result uses the multiplier method. Part of the difficulty is caused 
by the fact that, as shown in [l], there is no differential multiplier that provides us with 
a positive local energy norm. What we do instead is find a suitable pseudodifferential 
operator that does the job. This is chosen so that its symbol vanishes on trapped rays, 
which leads to a local energy norm which is degenerate there. 

As in the Schwarzschild case, the choice of re € (r_,r+) is unimportant since the 
r-slices r = const g {r-,r+) are spacelike. Hence in the proof of the theorem one can 
assume without any restriction in generality that is close to r+. 

Proof. The theorem is proved using a modification of the arguments in |21j . Let us first 
quickly recall the key steps in the proof of Theorem 13.11 as in [H]. We begin with the 
energy-momentum tensor 

Its contraction with respect to a vector field X is denoted by 

P4u,X]=Qal3[u]X'^ 

and its divergence is 

W"Pa[u,X] = DgU-Xu+^Qo^p[u]7T"f' 

where tt"'^ is the deformation tensor of X , given by 

TTq/S = ^aXp + V (}Xa 

A special role is played by the Killing vector field 

K = da 
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whose deformation tensor is zero. 

Integrating the above divergence relation for a suitable choice of X does not suffice in 
order to prove the local energy estimates, as in general the deformation tensor can only 
be made positive modulo a Lagrangian term of the form qd"udaU. Hence some lower 
order corrections are required. For a vector field X, a scalar function q and a 1-form m 
we define ^ ^ 

Pa [u, X, q, m] = Pa[u, X] + qudaU — —daqu^ + —rrLaU^. 
The divergence formula gives 

(4.12) V°'Pa[u,X,q,m] = Dgu(^Xu + qu'j + Q[u,X,q,m\, 

where 

Q[u,X,q,m] = ^Qap[u]Tr°''^ + qd°'udaU + maud°'u+{W°'ma - ^W°'daq)u^. 

So far these computations apply both for the Schwarzschild and the Kerr metrics. From 
here one, we will use the sub(super)scripts S, respectively K to indicate when a compu- 
tation is perform with respect to one metric or another. 

To prove the local energy decay in the Schwarzschild space-time, X, q and m are 
chosen as in the following lemma: 

Lemma 4.2. There exist a smooth vector field 

M 

X ^b{r){l^—)dr + c{r)K 
r 

with c supported near the event horizon and b > bounded so that 

Kb\ < Car-'' 

a smooth function q{r) with 

\d?q\ < Car-'-" 

and a smooth 1-form m supported near the event horizon r — 2M so that 
(i) The quadratic form Q^[u,X,q,m] is positive definite, 

3M^ ^ 



(4.13) Q^[u,X,q,m]>r-^\druf+ (^1 - — j {r-^\dvu\^ + r-'lful"^) + r-'u'^ . 

(a) X{2M) points toward the black hole, X{dr){2M) < 0, and {m,dr){2M) > 0. 

The local energy estimate is obtained by integrating the divergence relation (|4.12p 
with X + CK instead of X, where C is a large constant, on the domain 

£> = {0 < -Tj < Wo, r > re} 

with respect to the measure induced by the metric, dVs = r^drdvduj. This yields 

(4.14) J Q^[u,X,q,m]dVs = - j nsu{{X + CK)u + qu^dVs + BDR^[u] 
where BDR^[u] denotes the boundary terms 

/v=vo „ 
{di, P[u, X + CK, q, m])r^drduj - / {dr, P[u, X + CK, q, m])rldvduj 
v=0 J 
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Using the condition (ii) in the Lemma and Hardy type inequahties, it is shown in [21] 
that for large C and r^. close to 2M the boundary terms have the correct sign, 

(4.15) BDR'^[u]<ciE[u]{Y.-^)-C2{E[u]{v^) + E[u]{Y.+)), ci,C2>0 

Consequently, by applying the Cauchy-Schwartz inequality for the first term on the right 
of l|4.14p we obtain a slightly weaker form of the local energy estimate l|3.12p , namely 

(4.16) i?M(S+) + sn^E[u]{v) + \\u\\l^^, < i?M(S^) + \\f\\lEW'- 

V ^ ^ 

where the weaker norm LEWg and the stronger norm LEWg are defined by 



Mr 



1-^ 



2 



respectively 



Mr \ / 

These norms are equivalent with the stronger norms LEg, respectively LEg for r in a 
bounded set. On the other hand for large r the Schwarzschild space can be viewed as 
a small perturbation of the Minkovski space. Thus the transition from l|4.16p to p.l2p 
is achieved in [21] by cutting away a bounded region and then using a perturbation of 
a Minkowski space estimate. This part of the proof translates without any changes to 
the case of the Kerr space-time. Our goal in what follows will be to establish the Kerr 
counterpart of l|4.16p . namely 

(4.17) e[u]{j:+) + sup E[u]{d) + \\u\\l^,,n < e[u]{j:~^} + ll/llis^... 

where the norms LEW]^, respectively LEW^ coincide with LE]^, respectively LE'^ for 
bounded r, and with LEWg, respectively LEWg for large r. More precisely, if x('') is a 
smooth compactly supported cutoff function which equals 1 say for r < AM then we set 

IMlew^ = IMIIe]. + 11(1 - xMIew^ 

respectively 

Iklliw- = \\xu\\Ie'^ + - x)u\\Iew^ 
Different choices for x lead to different but equivalent norms. 

It is useful to first consider the effect of the same multiplier in the Kerr metric. The 
two metrics are close when measured in the same euclidean frame x = ruj with r > r^. 
Precisely, with d standing for dt and dx, x = ruj, 

(4.18) \d-[{gKh - (gsM < cc^, m{9Kf ~ {gsT']] < 
From this and the size and regularity properties of X, q and m it follows that 

(4.19) \P^[u,X,q,m]-P^[u,X,q,m]\ < ^\Vu\^ 
respectively 

(4.20) \Q^[u,X,q,m]-Q^[u,X,q,m]\<af^\Vu\^ + ^\u\ 
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Hence integrating the divergence relation (|4.12p in the Kerr space-time over the same 
domain D but with respect to the Kerr induced measure cIVk — p^drdvdio we obtain 

(4.21) j Q'^[u,X,q,'m\dVK = - J nKu({X + CK)u + qu^dVK + BDR'^[u] 

The bound l|4.19p shows that for small a the boundary terms retain their positivity 
properties in l|4.15p . namely 

(4.22) BDR''[u]<ciE[u]iE]^)-C2{E[u]iio)+E[u]{E+)), ci,C2>0 
However, (|4.20p merely shows that 

(4.23) QK[u,X,q,m]>r-^\dru\^ + 

and the right hand side is no longer positive definite near r = 3M. Thus we cannot close 
the argument as in the Schwarzschild case. As shown in [T], changing the vector field X 
near r ~ 3M would not help. 

To remedy this, we need to use a pseudodifferential modification S of the vector field 
X. We will choose S so that its kernel is supported in a small neighborhood of (3Af, 3M); 
this insures that there will be no additional contributions atr = r^. Furthermore, in order 
to be able to carry out the computations near the initial and final surfaces v — 0,va we 
take S" to be a first order differential operator with respect to v. Similarly, we modify the 
Lagrangian factor q using a pseudodifferential correction E, which is also a first order 
differential operator with respect to v. 

We also need to choose a quantization which is consistent with the Kerr measure. Here 
we have a few choices which have equivalent results. For our selection we use euclidean- 
like coordinates x = ujr. Given a real symbol s its eucHdean Weyl quantization is 
selfadjoint with respect to the euclidean measure dV = r^drdto. However, in our case we 
need to work with the Kerr induced measure dVx = p^drdto. Hence we slightly abuse 
notations and redefine the Weyl quantization as 

s'"' ■=- s"^ - 
p r 

If s is a real symbol, then (re)defined above is a selfadjoint operator in L^(dyR-). 

Another issue which does not affect our analysis but needs to be addressed is that 
we are using pseudodifferential operators in an exterior domain {r > r^} and some care 
must be given to what happens near r = r^. To keep things simple, in what follows 
all operators we work with are compactly supported in the sense that their kernels are 
supported away from re and infinity; even better, supported in a small neighborhood of 
3 A/. 

In what follows we consider a skewadjoint pseudodifferential operator S and a selfad- 
joint pseudodifferential operator E of the form0 

5 = < + s^9t, E = e]^ + -e_idt 



1 - 



3M 



Ca 



ir-^\diu\^+r-'\fu\^)+r- 



Since we are away from the event horizon the variable v coincides with t. We make this substitution 
here and later. 
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where si e , SQ,eo G 5*° and e_i e S^^ are real symbols, homogeneous outside a 
neighborhood of 0. Commuting and integrating by parts we obtain the counterpart of 
the relation l|4.2ip . namely 

(4.24) IQ'^[u,S,E] = I Oku- {S + E)udVK + BDR'^[u,S,E] 

Jd 

Here BDR^[u, S, E] represents the boundary terms obtained in the integration by parts 
with respect to t. Its exact form is not important, as all we need to use here is the bound 

(4.25) \BDR^[u,S,E] \ < E[u]{0) + E[u]{io) 

We note that due to the presence of the cutoff function x in both operators S and E 
there are no contributions on the lateral boundary r = r^. 

On the other hand IQ^[u, S, E] represents a quadratic form in (w, Dfu) which can be 
written in the form 

(4.26) IQ^^[u,S,E]^ f Q^u-u + 2^Q'^u-D^ + Q^Dtu'D^dVK 

Jd 

where QJ e OPS^ are selfadjoint pseudodifferential operators so that 

(4.27) + 2Qr A + Oo A' = \{PK, S] + U^E + EU^) 

We remark that for arbitrary operators S and E the expression above on the right is 
in general a third order differential operator in t. However, the operator S will always 
be chosen so that the coefficient of vanishes. We define the principal symbol of the 
quadratic form IQ^[u, S, E] as 

q''[S,E] = q2 + 2qiT + qoT^ 

The previous relation shows that it satisfies 

q^lS, E] = ^{p, s} + pe mod 5" + tS^^ + t^S'"^ 
2t 

We add (|4.2ip with a times (|4.24p . The boundary terms are estimated by l|4.22p and 
(|4.25p . Using the duality between the spaces CiL^ and L^. we can also estimate 

/ f{X + CK + q + a{S + E))udVK 
Jd 

Hence in order to prove l|4.17p it would suffice to show that the symbols s and e can be 
chosen so that 

(4.28) / Q''[u,X,q,m]dVK + aIQ''[u,S,E]>\\u\\l^^, 
Jd 

Here we aim to choose S and E uniformly with respect to small a. In effect, our construc- 
tion below yields symbols s and e which are analytic with respect to a. We remark that 
the choice of S and E is only important in the region where r is close to 3M. Outside this 
region, Q^[u,X,q,m] is already positive definite and the contribution of aIQ^[u^S^E\ 
is negligible. 

We consider first the expression Q^[u, X, q, m]. Near r — 3M this has the form 

(4.29) [u, X, q,m]=J2 q'^'^^dc^udpu + q'^^-^u^ 



< 11/11 



LEW' \\U\\ LEWI 
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where its principal symbol (7^^ = q^'°''^''laVi3 and the lower order coefRcient g^'" are given 
by the relation 

g^' = ^{p,X} + qp, <z^^" = -lv"d^q 

We do not need to exactly compute the above expression in the Kerr case, but it is useful 
to perform the computation in the simpler case of the Schwarzschild space. There we 
have 

where A stands for the spherical Fourier variable. Hence we obtain 

^4 r^^f - ^^r'v. ^} + (9 - r-H{r)(T - 3M)){r'p) 

where, near r = SM, 

rb{r){r - 3M)^ 



0i,{r) = ^b{r^ - 2Mr) ^) {b'{r^ - 2Mr) - b{r - A/)) 



(r - 2A/)2 

• ■ y r J 

respectively 

q{r) =q- r-'^b{r)ir - 3M). 

Here we have used the fact that 6 > to write the first two coefficients as squares. 

For our choice of q and r we know that the relation l|4.13p holds. This implies that 
the following two inequalities must hold: 

(4.31) g'5>e' + (r-3Af)2(T2 + A'), > 

Given the form of q^ , the first relation implies that g is a multiple of (r — 3Af)^, and that 
in addition there is a smooth function ^{r) so that 



r 



3 



2M 



q = iy{r)ag{r), < < 1 



This allows us to obtain the following sum of squares representation for q^ : 
(4.32) 

= (1 - iy{r))aliry + /32(r)f + i/i(r)a|(r)(A2 + (r^ ~ 2rA^)e'), i^i = v 

The symbol A^ of the spherical Laplacian ca also be written as sums of squares of differ- 
ential symbols, 

A = \-Y A2 ^~ A3 

where in Euclidean coordinates we can write 

{Ai, A2, A3} = {xi?7j - Xjr\^, i ^ j} 

(4.33) r^g^ = (1 - Hr))aUry + (iUr)e + v,{r)al{r){\\ + A^ + A^ + (r^ - 2rM)e) 
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We return now to the question of finding symbols s and e so that the bound (|4.28p 
holds. Near r ~ 3M, the principal symbol of the quadratic form on the left in l|4.28p is 

^{p,X + as} +p{q + ae) 

In order to prove (|4.28p at the very least we would like the above symbol to be nonneg- 
ative, and to satisfy the bound 

^{p, X + as)+p{q + ae) > cI{t - + cjir - + 

However, such a bound would not a-priori suffice since translating it to operator bounds 
would require using the Fefferman-Phong inequality, which does not hold in general for 
systems. Hence we prove a more precise result, and show that the symbols s and e can be 
chosen so that we have a favorable sum of squares representation for the above expression, 
which extends the sum of squares l|4.33p to a 7^ 0. 

Lemma 4.3. Let a he sufficiently small. Then there exist smooth homogeneous symbols 
s G S\^^, e G 'S'^Q^, also depending smoothly in a, so that for r close to iM we have 
sum of squares representation 

/I \ A 



(4.34) p^[-{p^X + as} + p{q , j- 

where fij G S'^o„; + 'T'S^iom satisfy the following properties: 

(i) The decomposition l|4.34p extends the decomposition (|4.33p in the sense that 

1 i 1 1 1 

(/xi,/i2,/i3,/i4,M5,M6) = ((1 - iy)^asT, f3s£,,i^i asXi,i^l asX2,i'i asX3,i^? ctsO 

mod a{Sl^^ + tS'^^^) 

and 

(a) The family of symbols {/ij}j=i,6 is elliptically equivalent with the family of symbols 
(c2(t — Ti), ci(t — T2), ^) in the sense that we have a representation of the form 

/ C2(t-Ti) 

11 = Mv, v= ci{t - T2) 

V e 

where the symbol valued matrix M G -^^^^'^('S'^om) ^C'S maximum rank 3 everywhere. 

Proof. Setting q = q — 2{ln p, X} respectively e = e — 2{ln p, s} we compute 

(^^{P, X + as} + {q + ae)p^ = Yi^P^^' X + as} + {q + ae){p'p) 

We first choose the symbol s so that the Poisson bracket {p^p, X + as} has the correct 
behavior on the characteristic set p = 0. Recall that the symbol oiX is ir^^b[r){r—iM)£^, 
where the vanishing coefficient at 3M corresponds exactly to the location of the trapped 
rays. Its natural counterpart in the Kerr space-time is the symbol 

s(r, r, ^, $) ~ ir^^b{r){r - roir, $))C. 
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This coincides with X in the Schwarzschild case a = 0, and it is well defined and smooth 
in a for r near 3M and |$| < 4|r|. In particular it is well defined in a neighborhood of 
the characteristic set p = 0, which is all we use in the sequel. 

We use (|4.7p to compute the Poisson bracket 

S} = - {p^p),r-'b{r)ir - ro(r, $)) + aP^p)^dr {r-^b{r){r - ro(T, $))) 

= 2r-^6(r)i?(r,T,$)A-2(r- ro(T,$)) 

+ [2Aar {r-^b{r){r - ro(T, $))) - 2(r - M)r-^h{r){r - ro(T, $))] 

Since ro(r, $) is the unique zero of i?(r, t, $) near ?' = 3M and is close to 3M, it follows 
that we can write 

(4.35) ^{p'p,s} = a'{r,T,^)T^{r~ro{T,^)f +P^{r,T,^)e on {p 0} 

where a,/? G <S'/|om ^-^^ positive symbols. We note that in the Schwarzschild case the 
symbols a and /3 are simply functions of r, see the first two terms in l|4.30p . 

Unfortunately s is not a polynomial in r, which limits its direct usefulness. To remedy 
that we first note that 

s-{iT-H{r){T-iM)0€aSl,^ 
Hence by (the simplest form of) the Malgrange preparation theorem we can write 

\s=[r- MI)b{r)^ + a(si(r, 9, 6, $) + so{r, ^ 9, 9, $)t) + a7(T, r, 0, 6, $)p 

with si S •S'^oni) '^0 S 'S'/Joni ^'^'^ 7 € ^hom- Then we define the desired symbol s by 

S = Sl + SqT 

The Poisson bracket \{p^Pi s} is a third degree polynomial in r. Hence, after division by 
p = —g**{T — ti){t — T2), taking also (|4.33p into account, we can write 

(4.36) tttIp^P) X + as} + Qip'^P) =12+ 1it + [es + a(eo + e_iT)](T - ti)(t - T2). 
where, by (j4.33p . the coefficient es corresponding to the Schwarzschild case is given by 

es = (1 - v{r))c?[r). 

It remains to show that the right hand side of l|4.36p can be expressed as a sum of 
squares as in the lemma modulo an error a(S^^^ + tSJ^^^Jp. Note that the symbols eg 
and e_i play no role in this, as they can be included in the error. 

The coefficients 71 and 72 can be computed using the relation l|4.35p and the fact that 
{p^p, X + as} = {p^p, 5} on p = (i.e. when t = Ti). We denote 

a, = J^a(r,T„$)(r-ro(T„$)) G 5^™, A = /3(r, r„ $), 

n - T2 

observing that ai can be used as substitutes for the c^'s in the lemma since they are 
elliptic multiples of Ci . Then we have 

l.yp,S}{n)='-aUn~T2)' + (3U' 
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which gives the following expressions for 71, 72: 
(4.37) 

72 = 7^1 - T2)(a2n - aiT-2) H $ , 71 = 7(^1 - T2)(ai - a^) + ^ 

We use the first components of 71 and 72 to obtain a sum of squares as follows: 
(ti - T2){alTi - a\T2) + t(ti - r2)(Q;i - al) = i^(ai(T - r2) - a2{T - Ti)f 
(4.38) + (1 - z/)(ai(T - T2) + a2(r - ti))^ 

- 4e/<-(T - ti)(t - T2) 

where 

(ai - az)^ , . 
eif = 1- (1 - y)aia2 

We remark that in the Schwarzschild case we have T2 = — ti and also ai = a2 = as and 
Pi ~ P2 ~ Ps- In particular this shows that 

eK- es € a(S'°o„ + tS^I^) 

which accounts for the eg factor in l|4.36p . It remains to consider the terms in (|4.37p . 
This is easier since the coefficients (32 are positive and have a small difference /3i — /32 G 
^^hom- Precisely, for a large C we can write 

ri/3|-T2/3? 1.^2^^2 ^ . , (Ca-/3| + /?i^)(r-r2)^ 
+ = 7;{Pi+ P2- Ca) + 

n - T2 n - T2 2 2(ti - r2)^ 

+ ^R^^^)^ + 

Summing this with l|4.38p we obtain the desired sums of squares representation, 

^{P^P,^ + as] + {p'^p)q e ^(ai(T - T2) - a2(r - ti))^ 

+ ^(«i(r - r2) + a2(r - n))' + + /^z' - Ca)e 

{Ca-l3l + l3l){T-r2)\ 2 , (Ca - + /3|)(r - n)^ , 
2(ri-r2)2 ^ ^ 2(ri-T2)2 ^ 

+ a(5Ln + ^w^)(T-Ti)(r-r2) 

Then e is chosen so that the last term accounts for the contribution of e. 

Part (ii) of the lemma directly follows. For part (i) we still need to specify which are 
the symbols ^j. Precisely, we set 

- - T2) + a2{T - ri))2, ^ 1 (^2 ^ ^2 _ ^^)^2 

2 _ (C«-/3|+/3f)(r-r2)% 2 2 _ [Ca - fij + fil){T - t,)\ ^ 

1{T,-T2f ^' ''^^ 2(n-r2)2 ^ 

Finally for ^3,4,5 and we set 

\2 

2 "^1,2,3 '^Z ^ \ / \\2 

M3.4,5 = ;^2 + (,2 _ 2,A/)^2 4(^1 - - «2(r - T,)) , 
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respectively 



2 _ (r^ - 2rM)e 



It is easy to see that for a = these symbols coincide with the coresponding Schwarzschild 
symbols. The proof of the lemma is concluded. 

□ 

In what follows we use the above lemma to prove the bound (j4.28p and conclude 
the proof of the theorem. We begin with symbols s and e as in the lemma. These 
are homogeneous symbols, and we can make them smooth by truncating away the low 
frequencies. They are only defined near r = 3M, therefore some spatial truncation is 
also necessary. Let x be a smooth cutoff function supported near 3M which equals 1 in 
a neighborhood of 3M, chosen so that we have a smooth partition of unity in r, 

At first we define the truncated operators 

S = X5"'X, E = xe^'x 
This choice would yield an expression 5, E] with a principal symbol 



ifrinAS, E]^x^ [ ^Ap^ s} + pe ) + ^xs{p, x} 



For these choices of S and E we consider the expression IQ^[u, S, E] which is given 
by l|4.26p with , and Qq as in l|4.27p . We observe that in general we can only say 
that the right hand side of l|4.27p is of the form 

^{pK, S] + UkE + EUk) = + 2Qr A + Q^Df + Q-^_^Dl Qf e OPS= 

However, its principal symbol q^^^^J^S^E] is at most a second order polynomial in r. 
Hence by the Weyl calculus we can write 

]^{[nK,S] + nKE+EUK)-{q^„nc[S.E]y" e OPS^'+OPS-^Dt+OPS-^Dl+OPS-^Dl 

In particular this shows that Q^i G OPS~^. To eliminate this term we slightly adjust 
our choice of E to 

E^Xe^X^Cu.Dt 
where the operator e^^^ is chosen so that 

9 ^aux ~^ ^aux9 = Q -1 

This is possible since the coefficient g** of in p is a scalar function which is nonzero 
near r = 3M. Also as defined e^„^ € OPS~^ and has kernel supported near r — 3M. 

Having insured that the term does not appear, we divide IQ^ [u, S, E] into two 
parts, 

IQ'^ [u, S, E] = IQ^„nc[n, S, E] + IQ!^^, [u, S, E] 
where the main component is given by 

(4.39) IQ^„,Ju,S,E]^ [ Q^^u-u + 2^Q'^^u-^+Q^^Dtu'D^dVK 

JD 
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with operators Q2.p-, Qi.p and Qo.p defined by 

while the remainder is given by a similar expression with operators Q^^; Qi a and QfC^ 
whose principal symbols are supported away from r = 3Af . More precisely, we have 

Qla + 2Qr,a A + OSia A' ^ (i^"^^' ^ + O^'^"' ^ + OP^"' A' 

Hence, using the fact that the LEW^ norm is nondegenerate outside an 0{a) neighbor- 
hood of 3M we can bound in an elliptic fashion 

(4.40) S,E]\< \\u\\l^^^^ + II A«|l^-i^^^ 

where the last term on the right represents the H^^ norm of Dj u in a compact region 
in r (precisely, a neighborhood of 3M). 

In order to conclude the proof of the theorem we turn our attention to the bound 
(|4.28p . which we seek to establish with S and E replaced with S, respectively E. We will 
show that 

(4.41) / Q''[u,X,q,m]dVK + aIQ^„„M,S,E]>\\u\\l^^, - 0(a)|| Aw||^-i 
We decompose the left hand side of l|4.4ip into an outer part and an inner part, 

where 

LHSW^^,= f xlQ''[u,X,q,m]dVk 
Jd 

L775{4311)„= / x^Q''[u,X,q,m]dVk + aIQ^„„,[u,S,E] 
Jd 

For the first part we use the pointwise positivity of away from 3M (see l|4.23p ) to 
conclude that 

(4.42) LHS^^^^> f xl{r-^\^u\^ + r-M^)dVK 

Jd 

The second part is a quadratic form which for convenience we fully recall here (see (|4.29p 
and l(439l) : 



Jd 

+ a [ Ql^pU-u + 2^Q'i;^u-^+Q^pDtu'D^dVK 
Jd 

where the coefficients g^'"'^, q^-^ and operators Q^^ satisfy 
respectively 

Qlp + '^QlpDt + QlpD^t - X (^{P, 4 + pe 
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We carefully observe that in the two parts of the expression for Li?5' l|4.41|) .j^ the cutoff 
function x appears in different places. In the first part, it is applied before the differen- 
tiation, while in the second part it is applied before the pseudodifferential operator. It 
does not make much sense to commute at this point. In the first part, we would produce 
lower order terms which may significantly alter q^'^ . In the second part, we would lose 
the compact support of the kernels for the operators Q™p- 

Since s and e are chosen as in Lemma 14. 3^ it follows that the principal symbol for 
Lg5' l|4.4iP ;,^ admits the sum of squares representation (|4.34p . We want to translate this 
into a sum of squares decomposition for LHS ^AAl^ .-^^. However, some care is required 
due to the different positions of the cutoff x, as explained above. The symbols ^fc = /Xfc(a) 
are in general of pseudodifferential type. However, part (i) of the Lemma guarantees that 
in the Schwarzschild case they are of differential type. Consequently, we write 

^lk{a) = Mfc(O) + ^lk[a) ~ ^fe(O) 

and use this decomposition to define the pseudodifferential operators 

Mk = XMfc(0)(x,i?) + {^lk{a) - Mfc(0))"'x 

Then using the Weyl calculus it follows that for Lg6' l|4.4ip ,,^ we have the representation 

JD f, JD 

where the remainder terms satisfy Vj G aS^~'^. What is important here is that the 
remainder is of size 0{a). This follows from our choice of the operators Mk, which 
guarantees that when a = the remainder is zero. 

Combining the last relation with l|4.42p we obtain the bound 



/ X^r-2|Vu|2 + r-4|«|2 + ^ \Mku\^dVK < LHS^ + + \\Dtu\ 

JD f. 



where the last two terms on the right account for the remainder terms involving the 
operators R^, which can be bounded using norms of u and Dtu in a compact region in 
r, away from r = and r = oo. 

It is easy to see that the above left hand side dominates IIwIIlbvKj^- For r away one 
uses only the first two terms. On the other hand for r close to 3M we use part (iii) of the 
Lemma, which guarantees that the symbols ci(t — T2), C2(t — n) and ^ can be recovered in 
an eUiptic fashion from the principal symbols fik of Mk- Thus l|4.4ip is proved. Together 
with l|4.40p this shows that 



u 



Iew^ ^ / Q''[u,X,q,m]dVK + aIQ''[u,S,E] + Oia)\\Dtu\\l- 

^ JD 



The final step in the proof of (|4.28p is to establish that the last error term above is 
neghgible. We can account for it in an eUiptic manner. Precisely, for any compactly 
supported selfadjoint operator Q g OPS~^ we can use in a Lagrangian term and 
integrate by parts (commute) to obtain 

25R I {g*')-^UKU-QhidVK^\\QDtu\\l,+0{ \\QDtu\\LA\u\\l,^^^ + \\u\\l,^^^ 

J D 

+ E[um + E[u][vo)) 
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which leads to the elliptic bound 

\\QDtu\\l.< Ml, +\\aj,u\\l^, + E[u]{Q) + E[u]{%) 

and further to 

\\Dtuf ^ <\\u\\l, +\\UKuf ^ + E[u]{Q) + E[u]{io) 

^comp comp ^comp 

Thus (|4.28p follows, and the proof of the theorem is concluded. □ 

Note that Theorem 14.11 tells us, in particular, that if we start with an initial data 
(uo,Mi) £ X L'^ then u{v) G is uniformly bounded for all -C > 0. A natural 
question to ask is if this is also true for higher iJ" norms. For n > 1 we define 

I Q I < n 

respectively 

|a| <n 

The higher order energies are similarly defined, 

\a\<n I I < n 

We then have the following 

Theorem 4.4. Let n be a positive integer and u satisfy ~ f with initial data 

(wo, 1*1 ) e X i?" on Sj^ and f G LE^*{Mb.)- Then 

i?"+iM(s]+) + supi?"+iM(5o) + ^ hoWl... + \\ui\\%, + ll/lli^.. 

v>0 ^ 

Proof. We remark that by trace regularity results we have 

|Q|<n-l 

Since the initial surface is space-like, we can use the equation to derive all higher v 
derivatives of u in terms of the Cauchy data (uo, ui) and /, 

Thus it suffices to prove that for uo > we have 

(4.43) E-+'[u]i^+) + E-+'[u]ivo) + Ml^.,. < E'^+'[u]{^n) + 11/11 W 

We will prove this for n — 2, and the proof for the other cases will follow in a similar 
manner by induction. 



Since is a Killing vector field, we have □A'(9ciu) = daf. Then by Theorem 14.11 we 
obtain 

(4.44) E[d,u]i^+) + E[d,u]ivo) + WOMIeI ^ E^[u]{E],) + \\f\\l^,. 
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In order to control the rest of the second order derivatives we take advantage of the 
equation, which takes the form 

(4.45) (.9™^,, + + L)u = f 

where L is a spatial partial differential operator of order 2. This is most useful in the 
region where is time-like. Given e > 0, this happens in the region of the form r > 2M+e 
provided that a is sufficiently small. The fact that (?s is time-like is equivalent to the 
ellipticity of the spatial part L of Oj^. From l|4.44p we obtain at v = vq 

The operator L on the left is elliptic in r > 2M -I- e, therefore by a standard elliptic 
estimate we obtain 

£;[Vu](S^„ n {r > 2M + e}) < E[u]ivo) + E[diiu]{vo) + WfWhi^^.^^) 

A similar elhptic analysis leads to the corresponding local energy bound, 

ll^"lliE]^(A<Hn{r>2M+2e}) ^ MlE}^(Mr} + \\(^'"'^\\lE]^{Mr) + \\f\\h{MR) 

+ E[Vu]{^Q n {r > 2M + e}) + £;[Vu](Eio n {r > 2M + e}) 

where the last two terms account for the output of integrations by parts in v. 

We are left to deal with the case r < 2M + 2e, where g^^ is small and simply using the 
equation l|4.45p does not suffice. Let C{r) be a smooth cutoff function such that C = 1 on 
[re, r_|. -I- 2e] and C = when r > r+ + 3e. Then we need bounds for the function w = Cu, 
which solves 

aKW = C.f + [aK,C]u:=9 

The commutator above is supported in the region {2M + 2e < r < 2M + 3e} where we 
already have good estimates for u. Recall that in the region {r < 2M + 3e} the LE]^ 
and LE^ norms are equivalent with the H^, respectively norm. Hence it remains to 
prove that for all functions w with support in {r < 2M + 3e} which solve U^w = g we 
have 

(4.46) i?[VH(S+)+i?[V^«](i}o) + ||Vi«||2,,(_^^) <i?2[^,](i]^) + ||g||2^,^_^^^ 

This is an estimate which is locahzed near the event horizon, and we will prove it taking 
advantage of the red shift effect. 

Since is Killing, this bound follows directly from Theorem l4.1l for the daw component 
of Vw, 

(4.47) E[d,w]{E+) + E[d,w](.vo) + \\d,w\\hiD) < E>]CS],) + Ml^^n) 
Consider now the angular derivatives of if, duiW. We know that 

[□5,9^] = 

since the Schwarzschild metric is spherically symmetric. Hence by (|4.18p it follows that 
[□jf , duj] is a second order operator whose coefficients have size 0{a). Hence by Theorem 
14.11 we obtain 

(4.48) E[d^w]i^+)+E[dMi^,o) + \\dM\HHD) < E^[w]i^^) + \\g\\HiiD)+a\\^ MlmiD) 
We still need to bound drW. For that we compute the commutator 

[□/f , dr]w — —{drg^^')drrW + Tw 
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where T stands for a second order operator with no terms. The key observation, which 
is equivalent to the red shift effect, is that drg"^^ > near r = 2M. Thus for X, C and q 
as in Lemma [321 we can write the equation for drW in the form 

(□k - j[{X + CK) + q])drw = drg + Tw 

with T as above and most importantly, a positive coefficient 7. Because of this the 
operator 

B = aK--f[{X + CK) + q] 
satisfies the same estimate in Theorem 14. II as Dg for functions supported near the event 
horizon. Indeed, the same proof goes through as in Theorem 13.11 Writing the integral 
identity l|4.2ip for w we see the contribution of 7 is negative, therefore we obtain the 
inequality 

J Q^[w,X,q,m]dVK < - J (drg + Tw)({X + CK)w + qu?jdVK + BDR^[w] 

By (|4.23p the left hand side is positive definite for r < 2M + 3e. Using Cauchy-Schwarz 
for the first term on the right and (|4.22p for the second we obtain 

Since T contains no second order r derivatives, this leads to 
(4.49) 

E[drW]i^+) + E[drW]ivo) + || 1 1 i (C) < E^[w]{^r) + MhHD) + "ll V^,5 w|| 

Then the desired bound l|4.46p follows by combining (|4.47p . I|4.48p and (|4.49p with ap- 
propriate coefficients. 

□ 

As an easy corollary, one obtains from Sobolev embeddings the pointwise boundedness 
result. 

Corollary 4.5. If u satisfies Oj^u = in Mr with initial data (mo,mi) G x in 
S^, then 
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